Conformal Invariance in (2+l)-Dimensional Stochastic Systems 



O 
(N 



I 

o 



(N 

0^ 
O 

ON 

O 



X 



L. Moriconi^ and M. Moriconi^ 

' Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
CP. 68528, 21945-970, Rio de Janeiro, RJ, Brazil 
^Instituto de Fisica, Universidade Federal Fluminense, 
Avenida Litordnea s/n, 24210-340, Niteroi, RJ, Brazil 

Stochastic partial differential equations can be used to model second order thermodynamical phase 
transitions, as well as a number of critical out-of-equilibrium phenomena. In (2+1) dimensions, many 
of these systems are conjectured (and some are indeed proved) to be described by conformal field 
theories. We advance, in the framework of the Martin-Siggia-Rose field theoretical formalism of 
stochastic dynamics, a general solution of the translation Ward identities, which yields a putative 
conformal energy-momentum tensor. Even though the computation of energy-momentum correlators 
is obstructed, in principle, by dimensional reduction issues, these are bypassed by the addition 
of replicated fields to the original (2-|-l)-dimensional model. The method is illustrated with an 
application to the Kardar-Parisi-Zhang (KPZ) model of surface growth. The consistency of the 
approach is checked by means of a straightforward perturbative analysis of the KPZ ultraviolet 
region, leading, as expected, to its c = 1 conformal fixed point. 

PACS numbers: 64.60.Ht, 64.60.De, 05.70.Jk 



Recent developments in two-dimensional criticality 
have extended the application of conformal field theory 
(CFT) methods [l|, to contexts far beyond the usual 
realm of second order phase transitions. Navier-Stokes 
turbulence 0, 01; random deposition models 043 1 spin 
glas ses [lol \VV\ and even chaotic quantum wavefunctions 
[1^ [Tsj provide actual examples where conformality has 
been convincingly established. The theory of stochastic 
Lowner equations (SLE's) [l^ - [i6l h as been a major tool 
in the study of these problems 17[. The essential strat- 
egy, which has strong geometrical appeal, is to identify 
level curves of order parameters or dynamical fields to the 
conformally invariant random paths described by SLE's. 
Nevertheless the clear success of numerical results, it is 
important to stress that they have been obtained in a 
frankly heuristic way, so that further elaborations are 
still in order. 

The understanding of critical two-dimensional out-of- 
equilibrium systems under the light of conformal meth- 
ods is still in its initial stages. The hope, supported by 
the above mentioned examples, is that conformal invari- 
ance will play a crucial role (as it does in thermodynam- 
ics) in extending the very concept of universality classes 
to the broad set of self-organized critical models (two- 
dimensional turbulence, for instance, turns out to be in 
the universality class of critical percolation Q). 

We put forward in this work an alternative, and essen- 
tially analytical, approach to the problem of conformal 
invariance in (2-f l)-dimensional dynamical systems. As 
an interesting testing ground, we apply the method to the 
Kardar-Parisi-Zhang (KPZ) model of kinetic roughening 
@j [l3 J ^ prototype of great interest in out-of-equilibrium 
statistical mechanics, due to its numerous connections to 
subjects like directed polymers in random media. Burg- 
ers hydrodynamics, flame fronts, domain wall dynamics, 
etc [li]. 

Let (p = (t>{x, t) be a general scalar field defined in 



(2 + l)-dimensional spacetime. Assume that correlation 
functions of (f> can be computed with the help of a path- 
integral measure £)0exp(—5' [(/>]), not necessarily real. We 
are interested to know if spatial fluctuations of (/) at a 
given time instant, say t = 0, have an underlying CFT 
description. Taking (f>^{x) = (f>{x,0), we introduce the 
reduced action Slcj)"] - the "0'*-model" (to be contrasted 
to the original (2+l)-dimensional "(/)-model") - which, 
up to normalization factors, is obtained from 



(1) 



where the symbol J stands for path-integration subject 
to the boundary condition ^(a;, 0) — (f>^{x). 

One may wonder if it is possible to define an energy- 
momentum tensor for the two-dimensional ^'^-model, as 
a way of probing its conformal structure from the com- 
putation of short-distance expansions [l|. Since it is in 
general not known how to evaluate it would seem 

hopeless to concentrate any effort on the computation of 
the corresponding energy-momentum tensor. However, 
as we show below, the energy-momentum route to con- 
formality is still viable in a pragmatical sense. Following 
Noether's theorem, let Tfj, with i,j = 1, 2, be the second 
order tensor related to the shift of the action S[(j)] induced 
by time independent infinitesimal translations ej (x) [l^ , 
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As Ti[, does not depend on the time variable, we call it the 



"projected energy-momentum tensor" . 
to show that 



It is not difficult 



mx) ^ expiSm) / D<l>eM-SmTtAx) (3) 
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does indeed solve the translation Ward identities of the 
model, viz., 

N 

p=i 

(4) 

We will, thus, pursue the idea that Eq. Q can be used to 
define the renormalized energy-momentum tensor associ- 
ated to the (/j^-model [2l|. It is not possible to compute 
T^j in general, since it depends, according to (P) and 
on exact path-integrations. Observe, however, that ([3]) 
leads to 



J D(f>' exp(-2S'| 



D(t>iD4>2S[ri - 02] exp(-5[0i] - S[(f>2]) 



xT^^{x,,X2)Tt{x[,x',) 



(5) 



The fields and 02 in ([5]) are replicas of 0, which are 
path-integrated under the constraint 0i(x, 0) = 02(x, 0). 
Eq. ([5]) would be exactly what one needs for the compu- 
tation of the energy-momentum correlator, if the expo- 
nential factor in its LHS were given by exp(— 5[0'']). 

It is clear that the critical surface is modified due to the 
replacement of exp(— S'[0'']) by exp(— 2S'[0*]) in the path 
integration measure. However, since the dimensionality, 
symmetries and the form of interactions are strictly the 
same in both of these models, it is likely that the critical 
surface can be restored to its original shape through cou- 
pling constant redefinitions. In the case of self-organized 
critical dynamics (as in the KPZ model) it is natural to 
conjecture stability of the critical behavior against such 
critical surface deformations. 

As far as we are interested in criticality, Eq. ([5]) 
is worth of attention due to the computability of its 
RHS, which avoids the use of the exact dimensionally 
reduced action S'[0'']. Taking w = xi + 1x2 and = 
Til — ^22 -l- 2zTi2, the holomorphic component of the pro- 
jected energy-momentum tensor derived in the 0-model, 
conformality of the 0*-modcl will be indicated by 
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where c is the central charge [l| and the above expec- 
tation value is computed from the path integration over 
fields 01 and 02 in ([5]). The 1/8 factor in ([5]) is due to 
the fact that the energy-momentum tensor associated to 
the model with action 2^[0^] is 2Tf^ . 

In order to understand how the method works in prac- 
tice, we focus our attention, now, on its application to a 
specific problem: the (2+l)-dimensional KPZ model. 



dth = ud^h+^idc.h)^ 



(7) 



which describes the evolution of the height h of randomly 
deposited atoms on a planar substrate. The system is 
characterized by the surface tension v, the non-linear pa- 
rameter A of lateral growth, and the strength of ran- 
dom deposition D. The gaussian random field rj is taken 
to have zero mean and two-point correlation function 
{rj{x, t)Tj{x', t)) =6'^{x- x')Sit - t'). 
The Cole-Hopf transformation, 



= eyi\){\h/2v) , 

maps Eq. ([7]) into 



2gi](f> , 



(8) 



(9) 



where g = 2\^D/v^. We can write down, using the re- 
sponse functional formalism [l^, IH, [2^ , the generating 
functional of correlation functions for the replicas 0i and 
02 of model ©, 



X exp I -S + i / dt(fx{jp4}p + jp4>p) 



(10) 



= d'xdt[t<j>p{dt(f>p~d'^p) + g^;^;] , (11) 



where S" = 5^ + S*^ with 



— 



I / d2xx(x)[02(x,O)-0i(f,O)] 



(12) 



Above, x{x) is the Lagrange multiplier field which en- 
sures, as in (O, the constraint 0i(£, 0) — 02(a;, 0). The 
system is confined, in principle, to a box of dimensions 
L X L, and time integration is restricted to the finite in- 
terval —A < t < 0, with eventually L, A — ^ cx) (the order 
of taking limits can be important here). 

Assuming that /i = at initial time {t — — A), the KPZ 
description is effectively a perturbation of the Edwards- 
Wilkinson (EW) linear model at small spacetime scales, 
where the surface height fiuctuates around h — f24j. 
We can perform, accordingly, the substitution — > 1 + 
in the action 5*^, to get 



9t0 = 9^0 + + V9V<t' 



(13) 



Perturbative renormalization group flows computed from 
either the stochastic Eqs. (O or (fT3|) agree perfectly well 
[H, [H [ii-IiSl, and show that g = is an ultraviolet 
fixed point in d = 2. Any finite g flows in the infrared to 
the (still barely understood) KPZ strong coupling regime 

As usual, the generating functional can be split into 
quadratic and interacting parts. Considering the small 
scale model (|13p . exact integrations over and yield 
the non-perturbed functional 



Z^xcxp i ^ I dtd^x 



dt'd^x' 



(14) 
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where Jp {x, t) is the q component of the doublet 



currents. We obtain 



jp{x,t) 
jp{x,t)-{-l)Px{x)6{t) 



(15) 



and Aqr{x,x' ,t^t') denotes the matrix elements of the 
operator 



a2 - 94 



~dt + 
dt + 92 -2ig 



(16) 



that is, 



Ai2{x,x',t,t') = A2l{x,x',t',t) 



1 ejt'-t ) 

2^ \t'~t\ 



■ exp 



(f-f ) 
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A22{x,x',t,t') = - — ln|f-f' I 

TT 

r-t-t' 



(f-f) 



n2 
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(17) 



The complete generating functional is written, then, as 
Z[{J^}] = 



= exp < / dtd X 



2g—. g—. 

SjpSjp 6j^6j 



(18) 



Applying time-independent spatial translations to the 
(/)-model's action, the holomorphic projected energy- 
momentum tensor is readily computed as 



T{w) = Ani J dtdu 



(19) 



The time variable is extended here to — 00 < i < 00. Note 
that due to the assumption of statistical stationarity one 
does not have to worry about boundary conditions at 
t -> ±00. 

Since the renormalized g vanishes at small scales, the 
exponential in (jlSp plays no role, and we expect the 
energy-momentum correlator to be given, in the ultra- 
violet region, by the free (quadratic) approximation: 

{Ti{w)T2{w')) = -{Anf j dt j dt' 



X [{dw(j)i{w, t)dw"j)2{w', t'))o{dw4>iiw, t)dnj'4>2{w', i'))o 
+ {du,$iiw,t)dw'(f>2iw' ,t'))o{d^(j)i{w,t)dw"p2iw' ,t'))o] . 

(20) 

The above two-point correlation functions are computed 
from the functional derivatives of Zo[{J|'}] at vanishing 
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: — J d^k- exp[-fc2(|t| + \t'\) +tk-{x- x')] , 

h{w,t)4>2{w',t')) = -e{-t)e{-t') X 

/ d^kk^cxp[~k^{\t\ + \t'\)+ik- {x-x')] , 

OT^ J 

h{w,t)Mw' ,t')) = 

f d^keM'k^{\t\ + \t'\)+^k-{x-x')] . 



87r2 



(21) 



Note that (0i(w, i)(/)2(w', i')) = {hiw,t)(j)i{w' ,t')). We 
are led, in this way, after a straightforward but lengthy 
computation, to 



{T,{w)T2{w')) 
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drrJ2{r) 
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8{w- w')" ' 



(22) 



where J2 (r) is the second order Bessel function of the first 
kind. It is remarkable that (|22p does not depend on the 
coupling constant g at all (for finite g). Comparing ([6]) 
and (122]), we find out that in the EW regime of the KPZ 
model, surface growing is conformal, with unit central 
charge, in agreement with the numerical SLE results @. 

At larger spacetime scales, the EW regime breaks 
down, as indicated by anomalous exponents in the 
Family- Vicsek (FV) finite size scaling relation for the 

roughness, W = y^{{h- hy) = L°'F{t/L'') [ll,^. In 
two dimensions, EW exponents are exactly a = and 
z = 2 while for the KPZ model, a ~ 0.38 and 

z ~ 1.62 [ii 113. 

Unfortunately, we are unable to access, by means of 
perturbation theory, the strong coupling (large scale) 
regime of the surface growing model The numerical 
analysis of Ref. points out that the KPZ universality 
class is conformal as well, with vanishing central charge. 
It is a challenging problem to derive it along the analyti- 
cal formalism presented here. It is worth mentioning that 
there are at least two non-perturbative approaches that 
could be useful, in principle, to address evaluations of the 
energy momentum correlator in the KPZ strong coupling 
regime: one is the "exact renormalization group" tech- 
nique [13]; the other is the 1/iV expansion developed in 
p. 

To summarize, we have discussed the issue of confor- 
mal invariance in (2-|-l)-dimensional stochastic systems 
within a systematic field-theoretical framework. A di- 
mensionally reduced energy-momentum tensor is defined 
in two dimensions from the original dynamical model, as 
a solution of translation Ward identities. Our main work- 
ing hypothesis is that critical surfaces are not spoiled 
by the introduction of replicated auxiliary fields, which 
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avoid computational complications associated with di- 
mensional reduction. The approach opens the way for 
analytical studies of conformality in a broad class of out- 
of-equilibrium models. A performance test has been car- 
ried out for the short-distance (weak-coupling) regime of 
the KPZ model of surface growing, where the method is 
noted to work with perfection. 
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